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We perform coupled-cluster and diffusion Monte Carlo calculations of the energies of circular 
quantum dots up to 20 electrons. The coupled-cluster calculations include triples corrections and a 
renormalized Coulomb interaction defined for a given number of low-lying oscillator shells. Using 
such a renormalized Coulomb interaction brings the coupled-cluster calculations with triples correla- 
tions in excellent agreement with the diffusion Monte Carlo calculations. This opens up perspectives 
for doing ab initio calculations for much larger systems of electrons. 
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I. INTRODUCTION 

Strongly confined electrons offer a wide variety of com- 
plex and subtle phenomena which pose severe challenges 
to existing many-body methods. Quantum dots in par- 
ticular, that is, electrons confined in semiconducting het- 
erostructures, exhibit, due to their small size, discrete 
quantum levels. The ground states of, for example, 
circular dots show similar shell structures and magic 
numbers as seen for atoms and nuclei. These struc- 
tures are particularly evident in measurements of the 
change in electrochemical potential due to the addition 
of one extra electron, = [i(N + 1) — /.J.(N). Here 
N is the number of electrons in the quantum dot, and 
fi(N) = E(N) — E(N — 1) is the electrochemical poten- 
tial of the system. Theoretical predictions of An and the 
excitation energy spectrum require accurate calculations 
of ground-state and of excited-state energies. 

The above-mentioned quantum mechanical levels can, 
in turn, be tuned by means of, for example, the appli- 
cation of various external fields. The spins of the elec- 
trons in quantum dots provide a natural basis for repre- 
senting so-called qubits 1 . The capability to manipulate 
and study such states is evidenced by several recent ex- 
periments (see, for example, Refs. dtl). Coupled quan- 
tum dots are particularly interesting since so-called two- 
qubit quantum gates can be realized by manipulating the 
exchange coupling which originates from the repulsive 
Coulumb interaction and the underlying Pauli principle. 
For such states, the exchange coupling splits singlet and 
triplet states, and depending on the shape of the con- 
fining potential and the applied magnetic field, one can 
allow for electrical or magnetic control of the exchange 
coupling. In particular, several recent experiments and 
theoretical investigations have analyzed the role of effec- 
tive spin-orbit interactions in quantum dots^~— and their 
influence on the exchange coupling. 

A proper theoretical understanding of the exchange 



coupling, correlation energies, ground state energies of 
quantum dots, the role of spin-orbit interactions and 
other properties of quantum dots as well, requires the 
development of appropriate and reliable theoretical few- 
and many-body methods. Furthermore, for quantum 
dots with more than two electrons and/or specific values 
of the external fields, this implies the development of few- 
and many-body methods where uncertainty quantifica- 
tions are provided. For most methods, this means provid- 
ing an estimate of the error due to the truncation made 
in the single-particle basis and the truncation made in 
limiting the number of possible excitations. For systems 
with more than three or four electrons, ab initio meth- 
ods that have been employed in studies of quantum dots 
are variational and diffusion Monte Carlo 8 ' 10 ! 11 , path in- 
tegral approaches^, large-scale diagonalization (full con- 
figuration interaction)^ 3 - -15 ' 17 , and to a very limited ex- 
tent coupled-cluster theory^—. Exact diagonalization 
studies are accurate for a very small number of electrons, 
but the number of basis functions needed to obtain a 
given accuracy and the computational cost grow very 
rapidly with electron number. In practice they have been 
used for up to eight electron a 13 ' 14 ' 17 , but the accuracy 
is very limited for all except N < 3 (see, for example, 
Refs. I15II21I ). Monte Carlo methods have been applied 
up to N = 24 elect ron o 10 ' 11 . Diffusion Monte Carlo, with 
statistical and systematic errors, provide, in principle, ex- 
act benchmark solutions to various properties of quantum 
dots. However, the computations start becoming rather 
time-consuming for larger systems. Hartree^, restricted 
Hartree-Fock, spin- and/or space-unrestricted Hartree- 
Fock^ 5 . - — and local spin-density, and current density 
functional methods^— give results that are satisfactory 
for a qualitative understanding of some systematic prop- 
erties. However, comparisons with exact results show 
discrepancies in the energies that are substantial on the 
scale of energy differences. 

Another many-body method with the potential of pro- 
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viding reliable error estimates and accurate results is 
coupled-cluster theory, with its various levels of trunca- 
tions. Coupled-cluster theory is the method of choice in 
quantum chemistry, atomic and molecular physic s 18 ' 32 , 
and has recently been applied with great success in nu- 
clear physics as well (see, for example, Refs. I38l - l4ll ). In 
nuclear physics, with our spherical basis codes, we ex- 
pect now to be able to perform ab initio calculations of 
nuclei up to 132 Sn with more than 20 major oscillator 
shells. The latter implies dimensionalities of more than 
10 100 basis Slater determinants, well beyond the reach 
of the full configuration interaction approach. Coupled- 
cluster theory offers a many-body formalism which allows 
for systematic expansions and error estimates in terms of 
truncations in the basis of single-particle states^ 2 -. The 
cost of the calculations scale gently with the number of 
particles and single-particle states, and we expect to be 
able to study quantum dots up to 50 electrons without a 
spherical symmetry. The main advantage of the coupled- 
cluster method over, say, full configuration approaches 
relies on the fact that it offers an attractive truncation 
scheme at a much lower computational cost. It preserves, 
at the same time, important features such as size exten- 
sivity. 

The aim of this work is to apply coupled-cluster the- 
ory with the inclusion of triples excitations through the 
highly accurate and efficient A-CCSD(T) approach 3 ^ 3 ^ 
for circular quantum dots up to N = 20 electrons, em- 
ploying different strengths of the applied magnetic field. 
The results from these calculations are compared in turn 
with, in principle, exact diffusion Monte Carlo calcula- 
tions. Moreover, this work introduces a technique widely 
applied in the nuclear many-body problem, namely that 
of a renormalized two-body Coulomb interaction. In- 
stead of using the free Coulomb interaction in an os- 
cillator basis, we diagonalize the two-electron problem 
exactly using a tailor-made basis in the centre-of-mass 
frame.— The obtained eigenvectors and eigenvalues are 
used, in turn, to obtain, via a similarity transformation, 
an effective interaction defined for the lowest 10 — 20 os- 
cillator shells. These shells define our effective Hilbert 
space where the coupled-cluster calculations are per- 
formed. This technique has been used with great success 
in the nuclear many-body problem, in particular since 
the strong repulsion at short interparticle distances of 
the nuclear interactions requires a renormalization of the 
short-range par t 43 ' 44 . With this renormalized Coulomb 
interaction and coupled-cluster calculations with triples 
excitations included through the A-CCSD(T) approach, 
we obtain results in close agreement with the diffusion 
Monte Carlo calculations. This opens up many interest- 
ing avenues for ab initio studies of quantum dots, in par- 
ticular for systems beyond the simple circular quantum 
dots. 

This article is organized as follows. Section II in- 
troduces (i) the Hamiltonian and interaction for circu- 
lar quantum dots, (ii) the basic ingredients for obtain- 
ing an effective interaction using a similarity-transformed 



Coulomb interaction, then (iii) a brief review of coupled- 
cluster theory and the A-CCSD(T) approach, and finally 
(iv) the corresponding details behind the diffusion Monte 
Carlo calculations. In Section III, we present our results, 
whereas Section IV is devoted to our conclusions and 
perspectives for future work. 

II. COUPLED-CLUSTER THEORY AND 
DIFFUSION MONTE CARLO 

In this section we present first our Hamiltonian in 
Subsection III Al thereafter we discuss how to obtain a 
renormalized two-body interaction in an effective Hilbert 
space. In Subsection III CI we present our coupled-cluster 
approach, and finally in Subsection lll Dl we briefly review 
our diffusion Monte Carlo approach. 

A. Physical systems and model Hamiltonian 

We will assume that our problem can be described en- 
tirely by a non-relativistic many-electron Hamiltonian H, 
resulting in the Schrodinger equation 

H\V) = E\V), (1) 

with | being the eigenstate and E the eigenvalue. The 
many-electron Hamiltonian is normally written in terms 
of a non-interacting part Hq and and interacting part V, 
namely 

N N 

H = H + V = ^2%i + ^2^, 

i—l i<j 

where Hq is the (one-body) Hamiltonian of the non- 
interacting system, and V denotes the (two-body) 
Coulomb interaction. In general, the Schrodinger equa- 
tion (fT]) cannot be solved exactly. 

We define the reference Slater determinant |$o) as the 
ground state of the non-interacting system by filling all 
the lowest-lying single-particle orbits. Since we will limit 
ourselves to systems with filled shells, this may be a good 
approximation, in particular if the single-particle field is 
the dominating contribution to the total energy. The 
non-interacting Schrodinger equation reads 

H m =eo|$o), (2) 

where 

N N 

Hp = hi = ^2 [*» + fi con(ri)] ■ 

i=l i=l 

The terms ti and v con (ri) are the kinetic-energy opera- 
tor and the confining potential (from an external applied 
potential field) of electron i, respectively. The vector 
represents the position in two dimensions of electron i. 
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Due to the identical and fermionic nature of electrons, 
the eigenstates of Eq. @ are Slater determinants, with 
the general form 

|$) = \ijk ...m) = ajaj4--aL|0), 

with a} being standard fermion creation operators (and a 
being annihilation operators). The single-particle eigen- 
states \i) — a\\0) and eigenenergies Si are given by the so- 
lutions of the one-particle Schrodinger equation governed 
by the operator hi. Since the total energy of the non- 
interacting system is given by the sum of single-particle 
energies e,-, we have 

N 

eo = X! £l ' 

i=i 

the reference determinant |<&o} is obviously the Slater 
determinant constructed from those orbitals with single- 
particle energies that yield the lowest total energy. In the 
particle-hole formalism, orbitals in the occupied space 
are referred to as hole states, while orbitals in the virtual 
space are denoted particle states. In principle, any com- 
plete and orthogonal single-particle basis can be used. 
However, since our coupled-cluster approach involves so- 
lution of a set of non-linear equations, it is preferable to 
start from a basis that produces a mean-field solution 
not too far away from the "exact" and fully correlated 
many-body solution. Therefore our main results will be 
obtained using the Hartree-Fock basis as a starting point 
for our coupled-cluster calculations. The Hartree-Fock 
basis is obtained from a linear expansion of harmonic os- 
cillator basis functions, such that the expectation value 
of the Hamiltonian is minimized. 

For the diffusion Monte Carlo calculations, it is also 
necessary to start from a model wavefunction that is used 
as importance function in the sampling, as we will dis- 
cuss later. The Slater determinant part, in this case, is 
built starting from the self-consistent orbitals generated 
in a Local Density Approximation calculation, in order 
to include as much information as possible about both 
exchange and correlation effects at the one-body level. 
Explicit two-body correlations are then included as an 
elaborate Jastrow factor; see Subsection III Dl for further 
details. 

Our model Hamiltonian^ for a quantum dot consists 
of a two-dimensional system of N electrons moving in 
the z = plane, confined by a parabolic lateral confining 
potential V con (r). The Hamiltonian is 

N fc2 2 N i 

" = »-^vf + Wr,)) + ^iir^ 

i— 1 i<j ' J ' 

In the above equation, m* is a parameter relating the 
bare electron mass to an effective mass, and e is the di- 
electric constant of the semiconductor. In the following 
(if not explicitly specified otherwise), we will use effec- 
tive atomic units, defined by h = e 2 /e = m e m* — 1. In 



this system of units, the length unit is the Bohr radius ao 
times e/m*, and the energy has units of Hartrees times 
m* /e 2 . As an example, for the GaAs quantum dots, typ- 
ical values are e = 12.4 and to* — 0.067. The effective 
Bohr radius a*, and effective Hartree H* are ~ 97.93 A 
and ~ 11.86 meV, respectively. In this work we will con- 
sider circular dots only with N = 2, N = 6, N = 12, 
and N = 20 electrons confined by a parabolic potential 
V con (r) = m e m*uj 2 r 2 /2. The numbers N = 2, N = 6, 
N = 12, and N = 20 are so-called magic numbers cor- 
responding to systems with closed harmonic oscillator 
shells, and hopefully a single- reference Slater determi- 
nant yields a good starting point for our calculations. 
The one-body part of our Hamiltonian becomes 

N / 1 2 \ 

whereas the interacting part is 



The unperturbed part of the Hamiltonian yields the 
single-particle energies 

e; = ui (2n + \m\ + 1) , (3) 

where n — 0, 1, 2, 3, .. and m = 0, ±1, ±2, ... The index 
i runs from 0,1,2,.... The shell-structure is clearly de- 
duced from this expression. We define R as the shell 
index. We will denote the shell with the lowest energy as 
R = 1, the shell with the second lowest energy as R = 2, 
and so forth. Hence, 

i?, = — (i= 1,2,3,...). (4) 

U! 

In the calculations, we limit ourselves to values of u) = 
0.28 a.u. (atomic units), uj — 0.5 a.u., and u) — 1.0 a.u. 
For higher values of the oscillator frequency, the contri- 
bution to the energy from the single-particle part dom- 
inates over the correlation part. The value uj — 1.0 is 
an intermediate case which also allows for comparison 
with Taut's exact solution for N = 2, see Ref. |46|, while 
uj = 0.5 and 0.28 represent cases where correlations are 
stronger, due to the lower average electron density in the 
dot. 



B. Effective interaction 

Whenever a single-particle basis is introduced in or- 
der to carry out a many-body calculation, it must be 
truncated. The harmonic oscillator basis is the de facto 
standard for quantum dots and nuclear structure calcula- 
tions. In nuclei, the intrinsic Hamiltonian is most easily 
treated using this basis, and for quantum dots the con- 
fining potential is to a good approximation harmonic. 
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However, the discrete Hilbert space R obtained from 
such a truncation grows exponentially with the number of 
particles. For example, allowing n single-particle states 
and N particles, 

dim(H) = (l\ . 

As an example, if we distribute N — 6 electrons in 
the total number of single-particle states defined by 10 
major oscillator shells, we have n = 110, resulting in 
dim('H) ~ 2.3 x 10 13 Slater determinants. This number 
is already beyond the limit of present full configuration 
interaction approaches. In our coupled-cluster calcula- 
tions we perform studies up to some 20 major shells. 
For 20 shells, the total number of single-particle states 
is n = 420, for which dim(R) « 1.3 x 10 18 , well beyond 
reach of standard diagonalization methods in the foresee- 
able future. For 20 electrons in 20 shells, the number of 
Slater determinants is much larger, 7.G x 10 ,w in total. 

But even if we could run large configuration interaction 
calculations, the convergence of the computed energies as 
a function of the chosen single-particle basis is slow for 
a harmonic oscillator basis, mainly due to the fact that 
this basis does not properly take into account the cusp 
condition at |r$ — Tj\ =0 of the Coulomb interaction. In 
fact, the error AE in the energy for a quantum dot prob- 
lem, when increasing the dimensionality to one further 
shell with a harmonic oscillator basis, behaves like 

AE-OiR^ 5 - 1 ). (5) 

Here, k is the number of times a given wave function 
\E' may be differentiated weakly, 8 G [0, 1) is a constant 
and Rho is the last oscillator shell. The derivation of the 
latter relation is detailed in Ref. |2l|, together with exten- 
sive discussions of the convergence properties of quantum 
dot systems. For the ground state of the two-electron 
quantum dot, we have precisely k — 1, while for higher 
electron numbers one observes k — O(l). This kind of 
estimate tells us that an approximation using only a few 
HO eigenfunctions necessarily will give an error depend- 
ing directly on the smoothness k. 

Although the coupled-cluster method allows for the 
inclusion of much larger single-particle spaces, the slow 
convergence of the energy seen in full configuration inter- 
action calculations applies to this method as well as it ap- 
proximates the configuration-interaction solution using 
the same set of single-particle functions. For an overview 
of coupled cluster error analysis, see Refs.HHl 

One way to circumvent the dimensionality problem is 
to introduce a renormalized Coulomb interaction V c r de- 
fined for a limited number of low-lying oscillator shells. 
Such techniques have been widely used in nuclear many- 
body problems (see, for example, Refs. [343.48). For 
quantum dots, this was first applied to a configuration 
interaction calculation by Navratil et aZ42, albeit for a 
different quantum dot model. But the potential of this 
method has not been explored fully, except for recent 



prelimary studies in Refs. [l5ll2lll50l . which demonstrate 
a significant improvement of the eigenvalues. Further- 
more, we expect that the potential of this method is of 
even greater interest when linked up with an efficient 
many-body method like the coupled-cluster approach. 

The recipe for obtaining such an effective interaction 
is detailed in several works (see, for example, Refs. Irsl 
[48ll56h . Here we give only a brief overview. 

The Hilbert space H is divided into two parts PH and 
QR, where P and is the orthogonal projector onto the 
smaller, effective model space, and Q = 1— P. Note here, 
that PR will be the space in which we do our many-body 
computations, and R is, in principle, the whole untrun- 
cated Hilbert space. The interaction operator V is con- 
sidered a perturbation, and we introduce a convenient 
complex parameter z and study H (z) = Hq + zV. Set- 
ting z = 1 recovers the original Hamiltonian. 

Consider a similarity transformation of H{z) defined 

by 

H{z) = e- x ^H{z)e x(z) (6) 

where the operator X(z) is such that the property 

QH(z)P = (7) 

holds. Equation ([6]) must not be confused with equations 
from coupled-cluster theory. The idea is that X(z) should 
be determined from perturbation theory, which gives an 
analytic operator function with X(0) = 0. 

The most important consequences of these equations 
are that (i) H have identical eigenvalues with H, (ii) that 
there are D = dim(P'H) eigenvalues whose eigenvetors 
are entirely in the model space PR. Thus, the effective 
Hamiltonian defined by 

H c g(z) EE PH{z)P (8) 

is a model-space operator with D exact eigenvalues. At 
z = these are the unperturbed eigenvalues, and these 
are continued analytically as z approaces z = 1. 

Equations ^ and ([7]) are not sufficient to determine 
X(z) uniquely. The order-by order-expansion of X{z) 
must be supplied with side conditions. One of the most 
popular conditions is that X(z)^ = — X(z) such that 
H(z) is Hermitian, and additionally that the effective 
eigenvectors are as close as possible to the exact eigen- 
vectors, i.e., that the quantity A defined by 

D 

A^^|||* fc )-|*e ff , fe )|| 2 (9) 

fe=l 

is minimized, where |\I/eff,fc) are the eigenvectors of -ff G ff, 
see Ref. [l5|. One can obtain a formula for X(z) in this 
case, namely 

X = artanh(w — a;'), 
where to = QujP is the operator such that exp(w)P|\l/fe) = 
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Order-by-order expansion of H(z) reveals that it con- 
tains ro-body terms for all m < N, even though V 
only contains two-body interactions. However, the many- 
body terms can be shown to be of lower order— 6 in z. By 
truncating H(z) at terms at the two-body level we obtain 
the so-called sub-cluster approximation to the effective 
Hamiltonian. This can be computed by exact diagonal- 
ization of the two-body problem; a simple task for the 
quantum dot problem^. 

The one-body part of H is always Hq 7 so it is natural 
to define the effective interaction by 

H cS = H + V cS . (10) 

The reader should, however, keep in mind that the sub- 
cluster approximation always produces missing many- 
body correlations when inserted in a many-body context. 
The size of this source of error can only be quantified a 
posteriori, either by comparison with experiment and/or 
exact calculations (see, for example, Ref. [H) for a dis- 
cussion on missing many-body physics and the nuclear 
many-body problem. 

C. Coupled-cluster method 

The single-reference coupled-cluster theory is based on 
the exponential ansatz for the ground-state wave function 
of the TV-electron system, 

|*o) = e T |$ ), 

where T is the cluster operator (an TV-particle- iV-hole ex- 
citation operator) and |$o) is the corresponding reference 
determinant (defining our chosen closed-shell system or 
vacuum) obtained by performing some mean-field calcu- 
lation or by simply filling the N lowest-energy single- 
electron states in two dimensions. 

The operator T is a simple many-body excitation op- 
erator, which in all standard coupled-cluster approxima- 
tions is truncated at a given (usually low) M-particle-M- 
hole excitation level M < N, with N being the number 
of electrons. If all excitations are included up to the N- 
particle- TV-hole set of Slater determinants, one ends up 
with solving the full problem. The general form of the 
truncated cluster operator, defining a standard single- 
reference coupled-cluster approximation characterized by 
the chosen excitation level M, is 

T = T 1 +T 2 + n + --- + T M , (11) 

where 

Tk = JkVfl E *ii."i°**oi ' ' •°L<* i * • • • <*»i ■ 
ii,...,t k ;ai,...,a* 

Here and in the following, the indices k, . . . label oc- 
cupied single-particle orbitals while a, b, c, . . . label unoc- 
cupied orbitals. The unknown amplitudes tf, tfj etc in 



the last equation are determined from the solution of the 
coupled-cluster equations discussed below. For a trun- 
cated T operator, we will use the notation T(M), where 
M refers to highest possible particle-hole excitations. 

As an example, we list here the expressions for 
one-particle-one-hole, two-particle-two-hole, and three- 
particle-three-hole operators, labeled T\, T 2 , and T3, re- 
spectively, 

T * = E E *? a & ( 12 ) 

i<Sf 0,>6f 

and 

T 2 = \ E E $ a ififoa*i (13) 

ij<Sf ab>€f 

and finally 

T * = ^ E E 1$gtfa&%&jk. (14) 
ijk<s / abc>sf 

We will in this work limit ourselves to a single reference 
Slater determinant $o- 

The cluster amplitudes t^"f n are determined by 
solving a coupled system of nonlinear and energy- 
independent algebraic equations of the form 

(4? 1 "f»|ff|$o) = > zi <•■•<;„, ai< •••<<!„, 

(15) 

where n = 1, . . . , M. Here, 

H = e- T ^He T ^ = (He T ^) c , (16) 

is the similarity-transformed Hamiltonian of the coupled- 
cluster theory truncated at M-particle-M-hole excita- 
tions and the subscript C denotes the connected part of 
the corresponding operator expression, and \^^"f n ) = 
a ai . . . a o-„ ai ^ . . . |$) are the n-particle-n-hole or n- 
tuply excited determinants relative to reference deter- 
minant |$o)- The Hamiltonians H and H are normal 
ordered. 

If we limit ourselves to include only one-particle-one- 
hole and two-particle-two-hole excitations, what is known 
as coupled cluster of singles and doubles (CCSD), the 
method corresponds to M = 2, and the cluster operator 
T( n > is approximated by 

T(M) = T(2) = J\ + T 2 , (17) 

given by the operators of Eqs. (fT!i|) and (flU]) . 

The standard CCSD equations for the singly and dou- 
bly excited cluster amplitudes t % a and t A, defining T\ and 
T2, respectively, can be written as 

($?|#(CCSD)|$ o )=0, (18) 

and 

<*#|JT(CCSD)|*)=0, i<j,a<b, (19) 
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where 

ff(CCSD) =H = e- T ^He T ^ = (He T ^) c (20) 

is the similarity-transformed Hamiltonian of the CCSD 
approach and the subscript C stands for connected dia- 
grams only. 

The system of coupled-cluster equations is obtained 
in the following way. We first insert the coupled-cluster 
wave function |\I>o) into the iV-body Schrodinger equa- 
tion, 

H\* ) = AE \* ), (21) 

where 

AE Q = E Q - {$ \H\<I> Q ) 

is the corresponding energy relative to the reference en- 
ergy ($o\H\$o), and prcmultiply both sides on the left 
by e~ T< ' to obtain the connected-cluster form of the 
Schrodinger equation 

i/|$)=A£ |$), (22) 

where 

H = e- T ^He T ^ = (H e T ^) c (23) 

is the similarity-transformed Hamiltonian. 

Next, we project Eq. (|22l) . in which T is replaced by 
its approximate form T(M), Eq. pip, onto the excited 
determinants l^" 1 '"""), corresponding to the M-particle- 
Af-hole excitations included in Tm- The excited determi- 
nants l^" 1 '""") are orthogonal to the reference determi- 
nant |$o)j so that we end up with nonlinear and energy- 
independent algebraic equations of the form of Eq. (IT51) . 

Once the system of equations, Eq. (|T5|) . is solved for 
t m or t*Y.'.'.o„ (or Eqs. (JUJ) and (JS) are solved for T a 

and T 2 or t % a and i^ b ), the ground-state coupled-cluster 
energy is calculated using the equation 

E = ($o|#|*o) + AE 

= (3 |#|*o) + ($o\B\$o}. (24) 

It can easily be shown that if H contains only up to 
two-body interactions and 2 < M < N, we can write 

E Q = ($ |# |$o) + (*o|[fl-(Ti + T 2 + ^Tf)]c\*o). (25) 



In other words, we only need T\ and T 2 clusters to calcu- 
late the ground-state energy E of the iV-body (N > 2) 
system even if we solve for other cluster components T n 
with n > 2. As long as the Hamiltonian contains up 
to two-body interactions, the above energy expression is 
correct even in the exact case, when the cluster opera- 
tor T is not truncated (see, for example, Refs. Il8ll32l for 
proof). 

The nonlinear character of the system of coupled- 
cluster equations of the form of Eq. fl~5|) does not mean 
that the resulting equations contain very high powers of 
Tm- For example, if the Hamiltonian H does not contain 
higher-than-pairwise interactions, the CCSD equations 
for the Ti and T 2 clusters, or for the amplitudes t l a and 
t l Q b that represent these clusters, become 

(*t\[H(l+T 1 +T 3 + ~I? + T 1 T 2 + ^lf)]c\*)=0, (26) 
<$#| [H(l + T X +T 2 + + T X T 2 + -7? 

J ZD 

+ \TlT 2 + ^Tf)] c |$) = 0. (27) 

The explicitly connected form of the coupled-cluster 
equations, such as Eqs. f(T5|) or (|26]l and (|27| . guarantees 
that the process of solving these equations leads to con- 
nected terms in cluster components of T and connected 
terms in the energy Eq, independent of the truncation 
scheme M used to define Tm- The absence of discon- 
nected terms in Tm and Eq is essential to obtain the 
rigorously size-extensive results^. It is easy to extend 
the above equations for the cluster amplitudes to include 
triples excitations, leading to the so-called CCSDT— hi- 
erachy of equations. Defining 

pq 

with f pq the Fock matrix elements and 

W = 2^2(pq\\rs){a+a+a r a s } 

pqrs 

where (pq\\rs) are anti-symmetrized two-body matrix el- 
ements, the extension to triples gives the following equa- 
tions for the amplitudes with one-particle-one-hole exci- 
tations 



(*i ILfTi + f(T 2 + 1/27?) + WT X + W(T 2 + 1/27?) + W{T X T 2 + 1/67? + T 3 )] c |*> = 0, 

I 



and with two-particle-two-hole excitations 
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($#| [/Ti + /(T 3 + T 2 Ti) + W + WT X + W(T 2 + 1/2T?) + W{T X T 2 + \/&Tf + T 3 ) 



+W{T 1 T 3 + 1/2T 2 2 + l/2T 2 T? + 1/242?)] |*) = 0. 

I 



and with three-particle-three-hole excitations we end up with 
I 

{^kllfn + f(T 3 Ti + 1/2T|) + WT 2 + W(T 3 + T X T 2 ) + W(1/2T 2 + T 3 T\l/2Tl + T{) 



+W(T 2 T 3 + l/2T 2 2 Ti + 1/2T 3 T 2 



1/6T 2 T?)] C |$> = 0. 



r 



Different approximations to the solution of the triples 
equations yield different CCSDT approximations. The 
CCSD method scales (in terms of the most computation- 
ally expensive contributon) as n^n^, where no represents 
the number of occupied orbitals and n u the number of un- 
occupied single-particle states. The full CCSDT scales as 

Coupled-cluster theory with inclusion of full triples 
CCSDT is usually considered to be too computation- 
ally expensive in most many-body systems of consid- 
erable size. Therefore triples corrections are usually 
taken into account perturbatively using the non-iterative 
CCSD(T) approach described in Ref. [H. Recently, a 
more sophisticated way of including the full triples known 
as the A-CCSD(T) approach, has been developed by 
Taube et alMrM.. In the A-CCSD(T) approach, the left- 
eigenvector solution of the CCSD similarity-transformed 
Hamiltonian is utilized in the calculation of a non- 
iterative triples correction to the coupled-cluster ground- 
state energy. The left eigenvalue problem is given by 

<*o|Afl- = JS?<* |A , (28) 

were A denotes the de-excitation cluster operator 

A = 1 + Ai + A 2 , 

with 

Ai = Ka a a\ , 



A2 = I H X ab a b a a4 a ] ■ 



(29) 

(30) 
(31) 



The unknowns, A^ and A^L result from the ground-state 



solution of the left eigenvalue problem (|28|) . Using a 
single-particle basis that diagonalizes the fock matrix / 
witin the hole-hole and particle-particle blocks simulta- 
neously, and utilizing the \ l a and A^ b de-excitation am- 
plitudes together with the cluster amplitudes, tf and tff , 



we get the non-iterative A-CCSD(T) energy correction 
to the coupled-cluster correlation energy (see Ref. for 
more details), 



AE 3 



1 



(3!) : 



]T (*o|A(/ hp + W)*H* 



ijkj 



ijkabc 



l 

abc 



x — m"k\(W N T 2 ) c \^) 



(32) 



Here, fh P denotes the part of the normal-ordered one- 
body Hamiltonian that annihilates particles and creates 
holes, while 



Tii'fe — fii fjj fkk faa fbb fa 



(33) 



is expressed in terms of the diagonal matrix elements of 
the normal-ordered one-body Hamiltonian /. In the case 
of Hartree-Fock orbitals, the one-body part of the Hamil- 
tonian is diagonal and fh P vanishes. The state I^Jfc) is a 
three-particle-three-hole excitation of the reference state. 
For a further discussion of various approximations to the 
triples correlations, see, for example, the recent review 
by Bartlett and Musiaii£. 

In this work we will focus on the CCSD, the CCSD(T), 
and the A-CCSD(T) approaches, using either a renormal- 
ized or an unrenormalized interaction. In order to avoid 
an iterative solution of the CCSD(T) and A-CCSD(T) 
equations, we start from a self-consistent Hartree-Fock 
basis such that the Fock matrix / is diagonal. Using such 
a basis, the computational cost of the CCSD(T) and A- 
CCSD(T) energy corrections is n^n^ number of cycles, 
done only once. It is also important to keep in mind, 
in particular when linking our coupled-cluster theory 
with Monte Carlo approaches, that a wavefunction-based 
method like coupled-cluster theory is defined within a 
specific subset of the full Hilbert space. In our case, the 
Hilbert space will be defined by all possible many-body 
wave functions which can be constructed within a certain 
number of the lowest-lying single-particle states. 
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D. Diffusion Monte Carlo 

The diffusion Monte Carlo method seeks the solution 
of the equation: 



t |*(R,t)> = [H - E ]\*{R,t)) 



(34) 



where R collectively indicates the degrees of freedom of 
the system (the 3N electron coordinates, in this case). By 
expanding the state ^(R, t)) on the basis of eigenstates 
\4> n ) of H, a formal solution of Eq. (|34p is given by 

|*(R,r)> - e-^- £ °) T |*(R,0)) 

= E„e-^-^) r |0„)(^|*(R,O)) (35) 

from which it is evident that for r —¥ oo the only surviv- 
ing component is the ground state of H . Eq. (j3"4"]l can be 
numerically solved by expanding the state to be evolved 
in eigenstates |Rj) of the position operator (called "walk- 
ers"), so that the evolution reads: 

^(R l |*(R,r))=^(R| e -^- B ^|R' i )(R' i |vl/(R',0)). 

i i 

Formally, in terms the Green's function solution of Eq. 
(l34l) . the solution can be written as: 



¥(R,t) = J G(R',R,r)*(R',0) dR'. 



(37) 



The Green's function G(R',R,r) = (R|exp[-(F- 
i?o)]|R) is in general unknown. However, in the limit 
At — > it can be written in the following form: 



G(R',R,r 




2nh 2 Ai 



(R— R') 

*7 



-[V(R)-£ ]At 



(38) 

that is as the product of the free particle Green's func- 
tion, having the effect of displacing the d-dimensional 
walkers, and a factor containing the potential, which is 
interpreted as a weight for the estimators computed at 
the walker position, and a probability for the walker it- 
self to generate one or more copies of itself in the next 
generation. Due to the divergence of the potential at the 
origin, it is necessary to modify the algorithm, introduc- 
ing the so-called "importance sampling" . In practice, the 
sampled distribution is modified by multiplying by an ap- 
proximate solution of the Schrddinger equation ^y(R), 
which is usually determined by a variational Monte Carlo 
calculation: 



* T (R)*(R,r) = J G(R', 



R. 



^*t(R')*(R',0) 

(39) 

A final important observation is the fact that the proce- 
dure described above is well defined only in the case of a 
totally symmetric ground state. For a many-Fermion sys- 
tem it would be necessary, in principle, to project on an 



excited state of the Hamiltonian, which leads to a severe 
instability of the variance on the energy estimation. This 
problem is usually treated by artificially imposing, as an 
artificial boundary condition, that the solution vanishes 
on the nodes of the trial function (fixed-node approx- 
imation) . Many other technical details enter the real cal- 
culation. A thorough description of the DMC algorithm, 
as implemented for the calculations of this paper, can be 
found in Ref. l52l . 

The fixed-node DMC calculations depend on the qual- 
ity of the trial wavefunction \&r(R), which is usually built 
starting from a parametrized ansatz. The values of the 
parameters are computed by minimizing the expectation 
value of the Hamiltonian on ^^(R). The trial wavefunc- 
tions we use have the formi^: 

Wconf 

*(R) i>g = exp[<t>{R)] Jl o,Sf > S (R) , (40) 



where the oti are variational parameters. Because in this 
paper we are considering only closed-shell dots which 
have L = and 5 = 0, the sum in Eq. (|4"D"j) reduces to a 
single term: 



s l=o,s=o = jy\ jyi 



(41) 



where the D x are Slater determinants of spin-up and 
spin-down electrons, using orbitals from a local density 
approximation calculation with the same confining po- 
tential and the same number of electrons. The function 
exp[0] in Eq. (l4Tfl) is a generalized Jastrow factor of the 
form: 



N 



m = E 



N i / 



.fc=i 



where 



l<3 



b(r) 



kirri 
aar. 



b(ri)r 



l + b(rj)rij 



(42) 



h 13 
°0 



b\ J tan" 1 [(r - R C ) 2 /2R C A] . (43) 



It explicitly includes one- and two-body correlations and 
effective many-body correlations via the spacial depen- 
dence of b(r). The quantity R c represents an "effective" 
radius of the dot, which is optimized in the variational 
procedure. The bo and b\ parameters depend only on 
the relative spin configuration of the pair ij. The pa- 
rameters Oj,- are fixed in order to satisfy the cusp con- 
ditions, that is, the condition of finiteness of the local 
energy H^/ty for r^- — > 0. For a two dimensional sys- 
tem, Oy = 1 if the electron pair ij has antiparallel spin, 
/and a,ij = 1/3 otherwise. The dependence of ay on the 
relative spin orientation of the electron pair introduces 
spin-contamination into the wavefunction. However, the 
magnitude of the spin contamination and its effect on the 
energy has been shown to be totally negligible in the case 
of well-optimized atomic wavefunctions^ and we expect 
that to be true here as well. 
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Also the coefficients jk in the one-body term, the co- 
efficients A, bo, and b\ in the two-body term, and the 
coefficients on multiplying the configuration state func- 
tions are optimized by minimizing the variance of the 
local energy^ 1 -. 

III. RESULTS 

We start our discussion with the results for the two- 
electron system, since these can, for certain values of the 
oscillator frequency, be compared with the exact results 
of Taut 46 . These results are presented in the next subsec- 
tion using both a renormalized two-body Coulomb inter- 
action and the 'bare' Coulomb interaction. Thereafter we 
present coupled-cluster results with singles and doubles 
excitations for systems with N = 6 and TV = 12 electrons 
with the bare Coulomb interaction. The slow conver- 
gence as a function of the number of oscillator shells with 
the bare interaction serves to motivate the introduction 
of an effective Coulomb interaction. In the main result 
section, we present CCSD, CCSD(T), and A-CCSD(T) 
results for N = 6,12 and N — 20 electrons using an ef- 
fective two-body Coulomb interaction and compare with 
diffusion Monte Carlo (DMC) calculations for the same 
systems. 



A. Results for two electrons 

In this subsection we limit our attention to the 
two-electron system and compare our DMC results 
with coupled-cluster calculations with CCSD correlations 
only. The results presented here serve to demonstrate the 
reliability of using an effective Coulomb interaction. 

The CCSD approach gives the exact eigenvalues for 
the two-particle system. We have employed a standard 
harmonic oscillator basis using the frequencies u> = 0.5 
and ui — 1.0 a.u. Our results are listed in Table fl] The 
variable R represents the number of oscillator shells in 
which the effective interaction case represents the model 
space for which the effective Coulomb interaction is de- 
fined. The calculations labeled CCSD-]/ represent the re- 
sults obtained with the unrenormalized or bare Coulomb 
interaction, while the shorthand CCSD-V^g- stands for 
the results obtained with an effective interaction. Since 
the latter, irrespective of size of model space (number 
of lowest-lying oscillator shells in our case) always gives 
the exact lowest-lying eigenvalues by construction (a sim- 
ilarity transformation preserves always the eigenvalues), 
these results are unchanged as a function of the num- 
ber of oscillator shells R. For the two-body problem, 
coupled-cluster theory at the level of singles and doubles 
excitations yields the same as exact diagonalization in 
the same two-particle space. In our case, the number of 
two-body configurations is given by all allowed configu- 
rations that can be constructed by placing two particles 
in the single-particle orbits defined by the given number 



of oscillator shells R. For lu — 1.0 a.u.,Taut's exact re- 





R 


CCSD-F 


CCSD-Ktf 


DMC 


0.5 


2 


1.786914 


1.659772 






4 


1.673874 


1.659772 






6 


1.667259 


1.659772 






8 


1.664808 


1.659772 






10 


1.660211 


1.659772 






12 


1.660091 


1.659772 






14 


1.660018 


1.659772 






16 


1.659970 


1.659772 


1.65975(2) 


1.0 


2 


3.152329 


3.000000 






4 


3.025232 


3.000000 






6 


3.013627 


3.000000 






8 


3.009237 


3.000000 






10 


3.000895 


3.000000 






12 


3.000654 


3.000000 






14 


3.000505 


3.000000 






16 


3.000406 


3.000000 


3.00000(3) 



TABLE I: Ground-state energies for two electrons in a cir- 
cular quantum dot within the CCSD approach with (CCSD- 
Ves) and without (CCSD-1/) an effective Coulomb interac- 
tion. The diffusion Monte Carlo (DMC) results are also in- 
cluded. For lu — 1 Taut's exact result from Ref. \^ is 3 a.u. 
All energies are in atomic units. There are no triples correc- 
tions for the two-body problem. The variable R represents 
the number of oscillator shells. 

suit from Ref. 13 is reproduced. The non-interacting part 
of the Hamiltonian gives a contribution of 2 a.u. to the 
ground-state energy while the two-particle interaction re- 
sults in 1 a.u. 

We notice also that the DMC results agree perfectly 
(within six leading digits) with our CCSD-14ff calcula- 
tions. The standard error in the DMC calculations is 
given in parentheses. 

If we, on the other hand, use the bare Coulomb interac- 
tion, we see that the convergence of the CCSD-V results 
as a function of R is much slower and in line with the 
analysis of Ref. |2l| and our discussion in subsection IIIBI 
One needs at least some 16-20 major oscillator shells (be- 
tween 272 and 420 single-particle states) in order to get a 
result within three to four leading digits close to the ex- 
act answer. The slow convergence of the bare interaction 
for the two-electron problem may be even more prevalent 
in a many-body system, in particular for small values of 
w, where correlations are expected to be more important. 
With more particles, we may expect even worse conver- 
gence. In Table |TT] we present for the case of u; = 1.0 
a.u. CCSD results for N — 6 and N — 12 electrons. 
The bare Coulomb interaction in an oscillator basis is 
used. The diffusion Monte Carlo results are for N = 6 
20.1597(2) a.u. and for N = 12 65.700(1) a.u. Using the 
bare interaction thus results in a slow convergence, as 
will be demonstrated in the next subsection. The result 
of 20.1742 a.u. obtained with an effective Coulomb at the 
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CCSD level for TV = 6 and R — 10 is much closer to the 
DMC result, as can be seen from Table IIIII These re- 



R 


N = 6 


N = 12 


2 


22.219813 




3 


21.419889 


73.765549 


4 


20.421325 


70.297531 


5 


20.319716 


66.989912 


6 


20.260893 


66.452006 


7 


20.236760 


65.971686 


8 


20.221750 


65.889324 


9 


20.211590 


65.838932 


10 


20.204345 


65.806539 



TABLE II: Ground-state energies for N = 6 and N = 12 elec- 
trons in a circular quantum dot within the CCSD approach 
using the bare Coulomb interaction. All energies are in atomic 
units. There are no triples corrections. Results are presented 
for an oscillator frequency cj = 1.0 a.u. The variable R rep- 
resents the number of oscillator shells. For N — 12 the first 
three shells are filled and there are no results for two shells 
only. 

suits serve the aim of motivating the introduction of an 
effective two-particle interaction. In the next subsection, 
we will make further comparisons between our results 
with and without an effective interaction. In particu- 
lar, we will try to extract convergence criteria for both 
approaches and link our numerical results with the pre- 
dictions made by Kvaal in Eq. ([5]). 

B. Results with an effective Coulomb interaction 

We present here our final and most optimal results for 
N = 6, N = 12, and N = 20 electrons using the CCSD, 
the CCSD(T), and the A-CCSD(T) approaches. We list 
the CCSD(T) triples results as well. This method has 
for a long time been considered as the calculational 'gold 
standard' in quantum chemistry due to its low computa- 
tional cost and accuracy. We emphasize, however, that 
the A-CCSD(T) approach is an improvement of the stan- 
dard CCSD(T) approach, and should therefore be consid- 
ered as our best and most accurate coupled-cluster cal- 
culation in this work. In all calculations we employ an ef- 
fective Coulomb interaction and a self-consistent Hartree- 
Fock basis for different values of the oscillator frequency 
lo and the model space R. The results are compared with 
diffusion Monte Carlo calculations (DMC)^. In addition 
to the values of u; = 1.0 and to = 0.5, which serve more 
as a reference for earlier calculations, we present results 
for to = 0.28 a.u., which corresponds to 3.32 eV, a fre- 
quency which should approximate the experimental sit- 
uation in Ref. H^. The role of correlations is also more 
important for smaller values of u, allowing us therefore 
to test the reliability of our single-reference CCSD and 
A-CCSD(T) calculations. As the system becomes more 
and more correlated, contributions from clusters beyond 



the T(3) (beyond three-particle-three-hole correlations) 
level might become non-neglible. For values of to > 1, 
the single-particle part of the Hamiltonian dominates and 
correlations play a less prominent role. 

Our results for N = 6, N = 12, and N = 20 electrons 
are displayed in Tables IIIII IIV1 and [V] respectively. We 
present also the mean- field energies (that is, the Hartree- 
Fock ground-state energies). These are labeled as Eq 
in the Tables. For all values of to with R — 20 ma- 
jor oscillator shells, our best coupled-cluster results, the 
A-CCSD(T) calculations, are very close to the diffusion 
Monte Carlo calculations. Even for 10 major shells, the 
results are close to the DMC calculations, suggesting 
thereby that the usage of an effective interaction provides 
a better starting point for many-body calculations. The 
convergence of the coupled-cluster calculation in terms of 
the number of major oscillator shells is also better than 
the results shown in Table [TTI with the bare Coulomb in- 
teraction. This discussion will be further elaborated at 
the end of this section. 





R 


E 


CCSD 


CCSD(T) 


A-CCSD(T) 


DMC 


0.28 


10 


7.9504 


7.6241 


7.6032 


7.6064 






12 


7.9632 


7.6245 


7.6023 


7.6057 






14 


7.9720 


7.6247 


7.6016 


7.6052 






16 


7.9785 


7.6249 


7.6012 


7.6048 






18 


7.9834 


7.6251 


7.6008 


7.6046 






20 


7.9872 


7.6252 


7.6006 


7.6044 


7.6001(1) 


0.5 


10 


12.1927 


11.8057 


11.7871 


11.7892 






12 


12.2073 


11.8055 


11.7858 


11.7880 






14 


12.2173 


11.8055 


11.7850 


11.7873 






16 


12.2246 


11.8055 


11.7845 


11.7868 






18 


12.2302 


11.8055 


11.7841 


11.7864 






20 


12.2346 


11.8055 


11.7837 


11.7862 


11.7888(2) 


1.0 


10 


20.6295 


20.1766 


20.1623 


20.1633 






12 


20.6461 


20.1753 


20.1602 


20.1612 






14 


20.6576 


20.1746 


20.1589 


20.1600 






16 


20.6659 


20.1742 


20.1580 


20.1592 






18 


20.6723 


20.1739 


20.1574 


20.1586 






20 


20.6773 


20.1737 


20.1570 


20.1582 


20.1597(2) 



TABLE III: Ground-state energies for N — 6 electrons in a 
circular quantum dot within various coupled-custer approx- 
imations utilizing an effective Coulomb interaction and the 
diffusion Monte Carlo (DMC) approach. The coupled-cluster 
results have been obtained with an effective two-body interac- 
tion using a self-consistent Hartree-Fock basis and the CCSD, 
the CCSD(T), and the A-CCSD(T) approaches discussed in 
the text. Eo is the Hartree-Fock energy while R stands for the 
number of major oscillator shells. All energies are in atomic 
units. 

In R = 20 major shells the A-CCSD(T) results are 
very close to the DMC results. As an example, con- 
sider the to = 1 results for N = 6 in Table HIH The 
CCSD result is 20.1737 a.u., while the A-CCSD(T) num- 
ber is 20.1582 a.u. The corresponding DMC energy is 
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20.1597(2) and very close to our A-CCSD(T) result. 
With R — 20 shells our coupled-cluster calculations are 
almost converged at the level of the fifth or sixth num- 
ber after the decimal point. At the end of this section we 
discuss the convergence properties of the various coupled- 
cluster approaches as functions of the number of oscilla- 
tor shells R. 

In Tables Hill El and [V] we see that the CCSD(T) re- 
sults are in most cases overshooting the diffusion Monte 
Carlo results. From numerous coupled-cluster studies in 
quantum chemistry, it has been found that CCSD(T) 
tends to overestimate the role of triples and thereby often 
overshoots the exact energy. The A-CCSD(T) approach 
has, on the other hand, been found to give highly accu- 
rate correlation energies, and even in some cases perform- 
ing better than the full CCSDT approach (see Refs.^^). 
This is also consistent with our findings for the CCSD(T) 
and A-CCSD(T) correlation energies in quantum dots. 

Let us briefly discuss the error in our coupled-cluster 
calculations. There are two sources of error, the first 
coming from the finite size of the single-particle basis, 
and the other from truncation of the cluster amplitude 
T at the T(3) excitation level (three-particle-three-hole 
exciations) . We are presently not able to provide a math- 
ematical error estimate on truncations in terms of the 
number of particle-hole excitation operators in the clus- 
ter operator T. However, several studies from quantum 
chemistry (see RefJ^ and references therein) and in nu- 
clear physic a 39 ' 41 have shown that the CCSD approach 
gives about 90% of the correlation energy while CCSDT 
gives about 99% of the full correlation energy. Assum- 
ing that the DMC results are to be considered as exact 
results, we can calculate the percentage of correlation en- 
ergy our CCSD and A-CCSD(T) calculations give for dif- 
ferent numbers of electrons N and values u> of the confin- 
ing harmonic oscillator potential. In Table [VTl we list the 
amount (in percentage) of correlation energy obtained at 
the CCSD and A-CCSD(T) level; the coupled-cluster cal- 
culations were done in a model space of R = 20 major 
oscillator shells. 

As we see from Table I VII the CCSD approximation 
gives 90%, or more, of the full correlation energy, while 
the A-CCSD(T) approximation is at the level of 99-100% 
of the full correlation energy for R = 20. The CCSD ap- 
proximation is clearly performing better for larger values 
lo of the confining potential, but this is expected since the 
system becomes less and less correlated for larger values 
of lo. This shows that our coupled-cluster calculations of 
circular quantum dots are within or even better than the 
accuracy seen in different applications in both quantum 
chemistry and nuclear physics. 

As previously discussed, DMC results reported in this 
paper are still affected by the fixed-node approximation. 
The extent of the error only depends on the nodal surface 
of the wavefunction. Because we use a single product 
of Slater determinants, given the circular symmetry of 
the dots considered, the nodes depend only on the set of 
single-particle functions used. Previous tests performed 



changing the set of single-particle orbitals show that dif- 
ferences are of the order of one millihartrees or less^. 
The optimization of the Jastrow factor only influences 
the variance of the energy, which is typically of the or- 
der of 0.5% of the total energy. Therefore, for circular 
quantum dots, we can conclude, assuming that the DMC 
calculations are as close as possible to the exact energies, 
that with an effective two-body interaction, a finite basis 
set of R = 20 major oscillator shells, and at most three- 
particle-three-hole correlations in the cluster amplitude, 
the remaining many-body effects are almost negligible as 
we are within 99-100% of the full correlation energy. 





R 


E a 


CCSD 


CCSD(T) 


A-CCSD(T) 


DMC 


0.28 


10 


26.3556 


25.7069 


25.6445 


25.6540 






12 


26.3950 


25.7066 


25.6388 


25.6491 






14 


26.4221 


25.7074 


25.6363 


25.6470 






16 


26.4410 


25.7081 


25.6346 


25.6456 






18 


26.4551 


25.7085 


25.6334 


25.6446 






20 


26.4659 


25.7089 


25.6324 


25.6439 


25.6356(1) 


0.5 


10 


39.9948 


39.2218 


39.1659 


39.1721 






12 


40.0409 


39.2203 


39.1599 


39.1667 






14 


40.0709 


39.2197 


39.1565 


39.1635 






16 


40.0922 


39.2195 


39.1543 


39.1615 






18 


40.1080 


39.2194 


39.1527 


39.1601 






20 


40.1202 


39.2194 


39.1516 


39.1591 


39.159(1) 


1.0 


10 


66.6596 


65.7552 


65.7118 


65.7149 






12 


66.7106 


65.7484 


65.7017 


65.7051 






14 


66.7445 


65.7449 


65.6961 


65.6996 






16 


66.7686 


65.7430 


65.6926 


65.6963 






18 


66.7867 


65.7417 


65.6903 


65.6941 






20 


66.8006 


65.7409 


65.6886 


65.6924 


65.700(1) 



TABLE IV: Same caption as in Table | 
are for A^ = 12 electrons. 



] except the results 



In order to study the role of correlations as a function 
of the oscillator frequency w and the number of electrons, 
we define the relative energy 



E-qmc — (Ho) 



E 



DMC 



(44) 



where (Ho) is the expectation value of the one-body op- 
erator, the so-called unperturbed part of the Hamilto- 
nian. For N = 6 this corresponds to an expectation 
value (Ho) = lOw for the one-body part of the Hamilto- 
nian, while for N — 12 and N — 20 the corresponding 
numbers are (Ho) = 28w and (Ho) = 60uj, respectively. 
Assuming that the diffusion Monte Carlo results are as 
close as possible to the true eigenvalues, the quantity e 
measures the role of the two-body interaction and corre- 
lations caused by this part of the Hamiltonian as func- 
tions of lo and N, the number of electrons. The results 
for e are shown in Fig. [T] Results for N — 2 are also 
included. 
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R 


E Q 


CCSD 


CCSD(T) 


A-CCSD(T) 


DMC 


0.28 


10 


63.2588 


62.2851 


62.1802 


62.1946 






12 


63.2016 


62.0772 


61.9503 


61.9692 






14 


63.2557 


62.0634 


61.9265 


61.9466 






16 


63.3032 


62.0646 


61.9214 


61.9423 






18 


63.3369 


62.0656 


61.9181 


61.9395 






20 


63.3621 


62.0664 


61.9156 


61.9375 


61.922(2) 


0.5 


10 


95.2872 


94.0870 


93.9864 


93.9971 






12 


95.3407 


93.9963 


93.8818 


93.8944 






14 


95.4164 


93.9921 


93.8700 


93.8833 






16 


95.4676 


93.9904 


93.8632 


93.8771 






18 


95.5043 


93.9895 


93.8588 


93.8730 






20 


95.5320 


93.9891 


93.8558 


93.8702 


93.867(3) 


1.0 


10 


157.4356 


156.0128 


155.9324 


155.9381 






12 


157.5613 


155.9868 


155.8978 


155.9042 






14 


157.6437 


155.9740 


155.8795 


155.8863 






16 


157.7002 


155.9669 


155.8687 


155.8758 






18 


157.7413 


155.9627 


155.8618 


155.8690 






20 


157.7725 


155.9601 


155.8571 


155.8646 


155.868(6) 



defined as 



TABLE V: Same caption as in Table II II I except the results 
are for N = 20 electrons. 







0.28 


L) - 


= 0.5 


U] = 


= 1.0 


N 


AE 2 


AE a 


AE 2 


AE a 


AE 2 


AE 3 


6 


94% 


99% 


96% 


100 % 


97% 


100% 


12 


91% 


99% 


94% 


100 % 


96% 


100% 


20 


90% 


99% 


93% 


100% 


95% 


100% 



TABLE VI: Percentage of correlation energy at the CCSD 
level (AE 2 ) and at the A-CCSD(T) level (AE 3 ), for different 
numbers of electrons N and values of the confining harmonic 
potential oj. All numbers are for R = 20. 



We see from this figure that the effect of the two- 
body interaction becomes increasingly important as we 
increase the number of particles. Moreover, the inter- 
action is more important for the smaller values of the 
oscillator frequency lo. This is expected since the con- 
tribution from the one-body operator is reduced due to 
smaller values of u>. Including more electrons obviously 
increases the contribution from the two-body interaction. 
Since our optimal coupled-cluster results are very close 
to the DMC results, almost identical results are obtained 
if we replace the DMC results with the A-CCSD(T) re- 
sults. 

We can also study the role of correlations beyond the 
Hartree-Fock energy Eq. In order to do this, we relate 
the Hartree-Fock energy Eq in Tables IHUfVl to the opti- 
mal coupled-cluster calculation, namely the A-CCSD(T) 
results. The relative difference between these quantities 
conveys thereby information about correlations beyond 
the mean-field approximation. This relative measure is 



X 



E 



A-CCSD(T) 



En 



E 



A-CCSD(T) 



(45) 



The results are shown in Fig. ^ior N — 6, N — 12, and 
N = 20. 

We see from this figure that correlations beyond the 
Hartree-Fock level are important for few particles and 
low values of uj. Increasing the number of electrons in the 
circular dot decreases the role of correlations beyond the 
mean-field approximation, a feature which can be under- 
stood from the fact that for larger systems, multi-particle 
excitations across the Fermi level decrease in importance. 
This is due to the fact that the single-particle wave func- 
tions for many states around the Fermi level have more 
than one node, resulting in normally smaller matrix ele- 
ments. Stated differently, with an increasing number of 
electrons, the particles close to the Fermi level are more 
apart from each other, in particular for those particles 
which occupy states around and above the Fermi level. 
The consequence of this is that correlations beyond the 
Hartree-Fock level decrease in importance when we add 
more and more particles. This means in turn that for 
larger systems, mean-field methods are rather good ap- 
proximations to systems of many interacting electrons in 
quantum dots. Similar features are seen in nuclei. For 
light nuclei, correlations beyond the mean field are very 
important for ground-state properties, whereas for heavy 
nuclei like 208 Pb mean-field approaches provide a very 
good starting point for studying several observables. 

The reader should, however, note that here we have 
limited our attention to ground-state energies only. 
Whether our conclusions about the role of correlations 
pertain to quantities like say spectroscopic factors re- 
mains to be studied. 

We conclude this section by studying in more detail the 
convergence properties of our coupled-cluster approaches, 
in particular, we will relate our A-CCSD(T) and CCSD 
results with the diffusion Monte Carlo results and study 
the dependence on R. This analysis will be performed 
with and without an effective Coulomb interaction. The 
reason for doing this is that we wish to study whether the 
convergence criterion of Eq. (|S|), derived for a full config- 
uration interaction analysis, applies to various coupled- 
cluster truncations as well. Furthermore, we wish to 
see whether our calculations with an effective interac- 
tion converge faster as a function of R compared to a 
calculation with the bare interaction. 

We compute the following quantities 



logi eccSD(-R) = log 10 
and 

l°SlO e A-CCSD(T) (R) = logio 



-^ccsd (R) — ^DMC 



E-* 



DMC 



(46) 



-Ea-ccsd(t) (R) - Eumc 



Ebmc 



(47) 
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In Fig. [3] we plot the results for N = 20 electrons and 
uj = 0.5. We have chosen these values since they repre- 
sent one of the cases where the A-CCSD(T) results are 
always above the DMC results and we have no crossing 
between these two sets of calculations. The CCSD re- 
sults on the other hand are always, for all cases reported 
here, above the DMC results. This means that the trend 
seen in Fig. [3] for the CCSD calculations applies to all 
cases listed in Tables |m]|V| while for the A-CCSD(T) cal- 
culations, these results are similar for all cases except for 
N = 6 and uj = 0.5 and uj = 1.0, N = 12, and oj = 1.0, 
and N = 20 and w = 1.0. In these cases the results at 
R = 20 are slightly below the DMC results. However, 
the agreement is still excellent. 

The interesting feature to note in Fig. [3] is that the 
CCSD results change marginally after R = 12 for N 
20. and there is essentially very little to gain beyond 
twenty major shells. With the present accuracy of the 
DMC results, we can conclude that the CCSD results 
reach at most a relative error of approximately 10~ 3 and 
that it stays almost stable from R = 12 shells. The 
relative error with respect to the Monte Carlo results 
does not change much. This applies to all CCSD results. 
This tells us clearly that there are important correla- 
tions beyond two-particle-two-hole excitations and that 
these correlations do not stabilize after some few shells. 
Furthermore, the slope of the A-CCSD(T) calculations 
is much more interesting and resembles the slope of the 
configuration interaction analysis of Ref. [2l| with an ef- 
fective interaction. For the ground states of three to five 
electrons, Kvaal found in Ref. [2ll a slope of approximately 
a = —4 to a — —5 for a parameterization 

log 10 e rj c + alog 10 i?, 

for the ground-state energies of various A^-electron quan- 
tum dots. The variable c is a constant. Our slopes vary 
between a = — 4 and a = — 6, resulting in a relative error 
of approximately 10" 5 at R = 20 for the results in Fig. [3] 
The slope of the A-CCSD(T) result is a = -4.93. The 
reader should note that the DMC results cannot reach 
a higher precision. The slope of the CCSD calculation 
with an effective interaction is a — —0.67 after R = 12. 

In the same figure, we plot also the A-CCSD(T) re- 
sults obtained without an effective Coulomb interaction, 
that is, with the bare interaction only. These results are 
labelled as log 10 £a-ccsd(t)(-R)— bare. A Hartree-Fock 
basis was used in this case as well in order to obtain con- 
verged solutions for the A-CCSD(T) equations. We see 
in this case that the convergence is much slower, result- 
ing in a slope given by a — —2.58, a result not far from 
the analysis of Ref. |2l| for the bare interaction. Figure [4] 
exhibits a similar trend, except that here we present re- 
sults for N — 12 electrons and uj = 0.5. The slope of the 
A-CCSD(T) results is now a = —6.38 with an effective 
interaction and a = —1.81 with a bare Coulomb inter- 
action. We notice again that the CCSD results saturate 
around R = 12 major shells. 

These results are very interesting as they show that 



the usage of an effective interaction can really speed up 
the convergence of the energy as a function of the num- 
ber of shells. Furthermore, these results tell us also that 
correlations beyond the singles and doubles approach are 
simply necessary. The convergence behavior of the A- 
CCSD(T) results resembles, to a large extent, those of a 
full configuration interaction approach with and without 
an effective interaction. Although we can extract similar 
convergence behaviors as those predicted in Ref. [2l| as 
functions of a truncation in the single-particle basis, the 
challenge is to provide more rigid mathematical conver- 
gence criteria for truncations in the number of particle- 
hole excitations. Here we can only justify a posteriori 
that triples corrections are necessary. Work along these 
lines is in progress. 



IV. CONCLUSIONS AND PERSPECTIVES 

We have shown in this work that coupled-cluster calcu- 
lations that employ an effective Coulomb interaction and 
a self-consistent Hartree-Fock single-particle basis repro- 
duce excellently diffusion Monte Carlo calculations, even 
for very low oscillator frequencies. This opens up many 
interesting perspectives, in particular since our coupled- 
cluster calculations are rather inexpensive from a high- 
performance computing standpoint. Properties like ad- 
dition spectra and excited states can be extracted us- 
ing equation-of-motion-based techniques (see, for exam- 
ple, Refs. Il8ll40ll4ll ) . Furthermore, since our codes run 
in an uncoupled basis, one can also study other trapping 
potentials than the standard harmonic oscillator poten- 
tial. A time-dependent formulation of coupled-cluster 
theory may even allow for studies of temporal properties 
of quantum dots such as the effect of a time-dependent 
perturbation. 

For circular dots, we found that with the inclusion of 
triples correlations, there are, for all systems studied, 
indications that many-body correlations beyond three- 
particle-three-hole excitations in the coupled-cluster am- 
plitude T, are negligible. We observe also that for 
systems with more particles, correlations beyond the 
Hartree-Fock level tend to decrease. Thus, although we 
are able to extend ab initio coupled-cluster calculations of 
quantum dots to systems up to 50 electrons, a mean-field 
description will probably convey most of the interesting 
physics. 

With two popular and reliable many-body techniques 
like coupled-cluster theory and diffusion Monte Carlo re- 
sulting in practically the same energies, one is in the posi- 
tion where one can extract almost exact density function- 
als for quantum dot systems. This allows for important 
comparisons with available density functionals for quan- 
tum dots. Finally, we have also noted that triples cor- 
relations are necessary in order to obtain correct results. 
The convergence pattern of our calculations resemble to 
a large extent those seen in full configuration interaction 
calculations. 
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FIG. 1: Relative correlation energy e defined in Eq. (144[) for 
different values of huj and number of electrons. The DMC 
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FIG. 2: Relative correlation energy \ defined in Eq. (|45[) for 
different values of hw and number of electrons. The numbers 
are obtained from Tables llllirvl using R = 20. 




FIG. 3: Relative correlation energy e defined in Eqs. (146 \ and 
(|47l) for different values of R. The values displayed here are 
for TV = 20 and cj = 0.5. The numbers are obtained from 
Table |V] We include also the A-CCSD(T) results obtained 
with the bare interaction. 
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FIG. 4: Relative correlation energy e denned in Eqs. (|46[1 and 
(|47l) for different values of R. The values displayed here are 
for TV = 12 and cj = 0.5. The numbers are obtained from 
Table IWl We include also the A-CCSD(T) results obtained 
with the bare interaction. 



